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Abstract: We investigate the relationship between super symmetric gauge theories with 
moduli spaces and matrix models. Particular attention is given to situations where the 
moduli space gets quantum corrected. These corrections are controlled by holomorphy. It 
is argued that these quantum deformations give rise to non-trivial relations for generalized 
resolvents that must hold in the associated matrix model. These relations allow to solve a 
sector of the associated matrix model in a similar way to a one-matrix model, by studying a 
curve that encodes the generalized resolvents. At the level of loop equations for the matrix 
model, the situations with a moduli space can sometimes be considered as a degeneration of 
an infinite set of linear equations, and the quantum moduli space encodes the consistency 
conditions for these equations to have a solution. 
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1. Introduction 

Recent developments on super symmetric gauge theories have shown that there is a deep 
relationship between SUSY theories that admit a large N limit and zero dimensional matrix 
models |i~7| , 19]. The connection between these seemingly different physical systems 



arose from the study of topological string theory amplitudes for open strings, but now 
there are purely field-theoretic arguments that show this connection. Of the two available 
proofs f?J of the correspondence (for a U(N) gauge theory), the most compelling one 
was given in the paper |5| , as it gave a full dictionary between loop operators on the matrix 
model and certain elements of the chiral ring of the gauge theory, as well as a proof that 
can be argued to be valid non-perturbatively. Part of the setup includes a correspondence 
between the gaugino condensate of the SUSY gauge theory and the 't Hooft coupling of 



the matrix model [19]. 



Having this connection makes it possible to solve for the vacuum structure of the 
supersymmetric gauge field theory if one knows how to solve the matrix model in the large 



N limit. We will refer the reader to [14, 31] for reviews of matrix models. A more recent 



list of solvable models appears in [38] 



The most successful and studied example to date has been given by deforming the 
U(N) M = 2 gauge theory, whose solution was given by Seiberg and Witten [^] to an 
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M = 1 gauge theory by adding a superpotential tr (V(X)), where X is the chiral field 
superpartner of the vector multiplet and V(X) is a polynomial in X, see for example 
g, [L9|, IH H H IH> H, 0, H- This theory is related to the one matrix model with 
potential The classical vacua of the gauge theory are described by distributing 

eigenvalues on the roots of the polynomial V'(X) giving rise to a theory whose classically 
unbroken gauge group is Y\ U(Ni), where there are iVj eigenvalues in the i-th. root of the 
polynomial V'(X). 

Similarly, in the matrix model, the vacua are also described by distributing (an infinite 
number of) eigenvalues over the classical saddle points of the matrix model. Each of 
these classical points in the eigenvalue plane becomes a cut in the eigenvalue plane due 
to quantum effects, when one considers the spectral curve of the matrix model. For each 
cut there is an associated 't Hooft coupling which counts how many eigenvalues end up 



in the given cut. This information is related to a partial gaugino condensate [19], whose 
holomorphic definition was given nonperturbatively in 0]. 

One of the advantages of this new method of computing gauge theory results is that 
it does not assume S-duality like symmetries to solve the theory, so it can provide tests of 
these highly non-trivial symmetries. For example, aspects of the N = 1* vacuum structure 
have been studied in this way p3| , |24{ ], but a general solution of the vacua of the theory in 
terms of matrix models is not yet available. 1 

Part of the failure to obtain the full solution of the matrix model comes from the fact 
that the classical M = 1* gauge theory vacua is described by a decomposition into irre- 
ducible representations of SU(2), and there are an infinite number of such representations. 
In the associated matrix model, this translates into an infinite number of classical saddle 
points for the eigenvalues of one of the matrices, let us call it X, and to each such classical 
saddle point for the eigenvalue x of X one can associate a cut in a spectral curve of the 
quantum theory (at least, this is what experience dictates from the one-matrix model). 
In particular, one has a potentially infinite number of cuts. This suggests that one can 
not solve the loop equations of the matrix model in a straightforward fashion, as there 
will be an infinite number of unknown parameters: one can take these parameters to be 
the number of eigenvalues in each cut Ni, but more to the point, the one point functions 
< X n > depend classically on all of the N{. Thus, there are no recursion relations which 
solve for the one point functions of the loop operators < X n >. 

The technical issue is how to effectively truncate the problem to finitely many cuts. 
In the one matrix model this is automatic, but in the M = 1* theory one has to do 
this by hand. A related problem in complex analysis can be phrased as follows: given a 
Laurent series of an analytic function f(z) = Y^ a i z ~ l at infinity, what are the constraints 
on the variables dj that guarantee that f(z) has a finite number of cuts in the interior of 
the complex plane? The constraints involve an infinite number of the coefficients of / at 
the same time, since multiplying f{z) by any polynomial function of z will always give a 



1 This theory can also be solved by first turning on a mass term deformation that gives rise to an Af — 2 
gauge theory, and finding the associated integrable system p2| , pt| . In principle, it is possible to evaluate 
all of the holomorphic information by going to the appropriate point in the moduli space of vacua. 



- 2 - 



function with the same properties. Looking at it from this point of view it also looks like 
the problem has no solution. 

The current literature solves the abstract problem by making a clever anzatz for the 
one-cut solution (or multi-cut solutions where all of the cuts are correlated, see for example 
1 39, 23, 24, |25], f45|), that leads to studying functions on an elliptic curve, but it does not 
solve in any sense a general multi-cut solution. 

In this paper we will describe a new technique to solve (only partially) similar matrix 
models based on getting information from the gauge theory first. Since the gauge theory 
and the matrix model data are in some sense equivalent, this statement might seem para- 
doxical. We will concern ourselves with situations where the superpotential is non-generic 
and there is the possibility of having a moduli space of vacua. In particular, the allowed 
deformations of the moduli space give a finite number of parameters which determine part 
of the vacuum structure. In the matrix model these turn out to be situations where the 
system of loop equations is degenerate in some sense, a statement which will be explained 
in detail later in the paper, and there are more free parameters than at a generic case. The 
main question we want to address is: what does the moduli space tell us about the matrix 
model? 

To avoid undue suspense, we find that under some circumstances the existence of a 
moduli space of vacua produces an integrable sector of the matrix model, this is, part of 
the problem behaves like a one matrix model: one can find an infinite number of one 
point functions with a finite amount of data. The finite amount of data is exactly the 
one that characterizes the possible deformations of the moduli space that are allowed by 
holomorphy. 

To show the new technique, we will study two examples. First, a toy example which 
reduces to the one matrix model, which we understand very well. It will just be an illus- 
tration of the matrix model techniques used. Then we will study a three matrix model 
closely related to the M = 1* gauge theory which has special properties that insure that 
certain objects in the matrix model only have a finite number of cuts automatically. This 
happy coincidence is tied to the fact that these particularly gauge theories can be geomet- 
rically engineered. We will in this paper always consider a gauge field theory which can 
be geometrically engineered by placing branes on a singular Calabi-Yau (CY) geometry. 
From this point of view we are exploring questions that relate the CY geometry to a matrix 
model. 

The basic idea is that if we place fractional branes at the singularities of the geometry, 
then the CY geometry will be deformed by geometric transitions so that the singularities 
with fractional branes on them are 'resolved'. This will be exactly like in the Klebanov- 
Strassler geometry |33|, repeated many times at many singularities, as in |l2|, [10[] . 

If we insert a probe brane in the bulk (in the presence of the fractional branes) then it's 
moduli space should correspond to the deformed CY geometry. This idea is very common 



in the literature, see for example |40|, |34j, |35|, |4J . Fortunately, the shape of the CY can be 



studied directly from the field theory || [| so it is possible to start with a classical gauge 
field theory and produce certain CY geometries, as opposed to start with the geometry and 
guess the field theory. If we place a brane in the bulk and compute it's moduli space, the 



- 3 - 



moduli space of the brane will be the CY geometry. Of course, not every superpotential 
will lead to a CY geometry, but for those that do one can make progress in understanding 
the gauge theory questions by studying the geometry (moduli space of the theory). Also, 
one should be able to understand aspects of the matrix model by studying the geometry. 

The main advantage with this setup is that the problem is controlled by holomorphy: 
the possible deformations of the CY geometry that result from placing fractional branes 
at singularities are not arbitrary, so there are only a finite number of unknowns. This 
is exactly what we are looking for, a situation where the potentially infinite number of 
variables (cuts in the matrix model) is reduced to finitely many (we are assuming that 
there is one cut per quantum deformation parameter, a statement which will be justified 
later) . 

We want to study exactly what matrix model information we can get from the gauge 
theory moduli space. We will use a proposal first presented in j| to calculate the quantum 
moduli space from the matrix model. The relation between the moduli space and how to 
obtain it from the matrix model still lacks a proof, so it should be treated as a conjecture 
that needs justification. However, one can make predictions based on the calculation of the 
deformed moduli space from the matrix model, and the one expected from holomorphy. 
Some of these predictions can be verified directly by manipulating the loop equations of 
the matrix model, and we will use this calculation as a consistency check of the proposal. 

In our examples, we will not solve the matrix model completely (find all one point 
functions from a finite amount of data). We will use the word 'solve' in the following more 
restricted sense: it is possible to determine an infinite number of one point functions of the 
matrix model with a finite amount of data. In essence, we obtain similar results to a one 
matrix model. 

The main example we will study is a gauge theory that will have three matrices X, Y, Z 
in the adjoint of U(N), and a superpotential of the general form 

tr (XYZ - qXZY + V(X)) 

with the restriction q n = 1, and a restriction on the form of V{X) which depends on n. 
Similar systems have been studied in ||, |?], |, |{J || || and are interesting on their 



own right as they encode structure of interesting deformations of M = 4 gauge theories 
1 44]. The restrictions placed on the superpotential guarantee that there is a CY geometry 



associated to the gauge theory [26, ^7], ||, and similar examples studied in |T^, O]. Also, 
after a linear change of variables on X, we can get for q ^ 1 an equivalent superpotential 
of the form 

tr (XYZ - qXZY + V(X) + m 2 Y Z) 
The M = 1* gauge theory has the same matter content and the potential 



tr 



(XYZ - XZY + m 2 X 2 + m 2 Y Z) 



so by varying q we can in principle get arbitrarily close to the M = 1* gauge theory. This 
particular model has been studied extensively, also in relation to AdS/CFT |l]]. This is 
the main physical motivation to write this paper: to approach the problem of solving the 
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J\f = 1* vacuum structure from the matrix models, with an arbitrary number of cuts, as well 
as the related q deformed systems. It is important to notice that setting q ^ 1 means that 
the conformal field theory with V = has only N = 1 SUSY in four dimensions. Therefore, 
techniques that depend on the quantum corrections being calculable in N = 2 theories, 
which then are softly broken to M = 1 can break down because there is less supersymmetry 
protecting the system. In particular, to use integrable systems like in Af = 2 gauge theories, 
but with less supersymmetry p6|] , one actually needs a proof that the moduli space does 
not get lifted by quantum corrections. These matrix models are also interesting on their 



own right, see for example [ 43 1 , so we can also gain insight into the problem of solving 



matrix models by using four dimensional SUSY gauge theory information. 
2. The one matrix model 

In this short section we will show how the analysis of a quantum moduli space can solve the 
one matrix model. The whole purpose of this section is to show that a quantum deformed 
CY geometry contains information about loop equations of the matrix model. 
Consider the CY threefold geometry 

uv = ( w -P{z)){w + P{z)) (2.1) 

where P(z) is a polynomial of degree n {P{z) = Y17=o aizt )> an( ^ which results from de- 
forming the A\ singularity geometry uv = w 2 . This procedure removes the codimension 
two singularities and leaves only codimension three singularities behind, which for generic 
P(z) are conifold singularities. 

The gauge theory associated to D-branes on this geometry has been studied in [12, 11, 
1C, ||, and it is given by a quiver diagram with two nodes, representing the afline dynkin 



diagram A\. 




Because there are two nodes in the quiver theory, one has general gauge group U (N) x 
U(M). The superpotential of the theory is given by 

W = tr (ZAiBx - B X A X Z - ZA 2 B 2 + B 2 A 2 Z) + tr (V(Z) - V(Z) (2.2) 

and in the case V{x) = the superpotential is the one that results from an M = 2 su- 
persymmetric gauge theory. P{x) as a polynomial is determined by the following equation 
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2P{x) = V'(x), so one can relate the shape of the geometry and the gauge theory param- 
eters in the superpotential. 

A brane in the bulk will have a theory characterized by the gauge group U(l) x U(l). 
This is the correct fractional brane content for a bulk brane in the orbifold geometry, and 
these numbers do not change when we deform the geometry. 

Also, if we take M = 0, we see that we obtain a gauge theory with the same matter 
content as the N = 2, U (N) gauge theory, which has been deformed by the potential V(Z). 
Thus this particular theory is a sub-theory of the one defined by the geometry. 

This U(N) does not have a moduli space of vacua, instead it has a discrete set of 
vacua. These vacua at the classical level are built by distributing the N eigenvalues of Z 
into the roots of P{Z). If there are N% eigenvalues at the rt root, the low energy effective 
(classical) theory is given by pure gauge theory, with gauge group \\iU{Ni). Quantum 
effects lead to confinement of the SU(Ni) groups, so in the infrared of the quantum theory 
we end up with a gauge group n(^(l)) over the roots which have at least one eigenvalue. 

To solve for the structure of the vacuum, we want to find the vevs of all of the elements 
of the chiral ring. The chiral ring is generated by the following combinations 

' tr{Z n W a W a ) (2.3) 



327T 2 

tr(Z n W a ) (2.4) 
tr{Z n ) (2.5) 

The ordering of W inside the trace does not matter, because any other ordering can be 
obtained by the addition of terms of the form DM, so at the level of the chiral ring they 
give the same cohomology class. For our purposes, the elements of the chiral ring which 
have a simple interpretation in a matrix model are the ones given by tr(Z n W a W a ). Using 
the generalized Konishi anomaly for the variations 5Z ~ Z k W a W a , Cachazo et al. (?J were 
able to show that these satisfy the same equations as the loop equations for the one matrix 
model with potential V(X). This is 

tr(V'(Z)Z k W a W a ) = — ^ ^(triZ^^WaW^triZ^W ) (2.6) 

32it z ^-^ 

If we call Oi = ^^jtr(Z l W a W a ), then the equations above read 

tr(J2 ajZ k+j W a W a ) = ^3 ^{triZ^-'WaW^triZiWaW ) (2.7) 
or equivalently 

a 3 O k+J = Ok^Oi (2.8) 

i 

This equation can be read as a recursion relation for the vevs O n +fc for k > in terms of 
the vevs Oq, . . . , O n -i- All of these equations can be put in one single equation in terms 
of a generating function for the Oj, defined as follows 

R(u) = Y J O k u- k - 1 (2.9) 
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so that 

V'{u)R{u) = R{uf + /(it) (2.10) 

where / is an unknown polynomial of degree n — 1. This is necessary to seed the initial 
conditions for the recursion relation. Given / it straightforward to solve for R{u). f in the 
quantum theory is determined by how we choose to distribute the eigenvalues of Z on the 
different roots of the classical potential, and it encodes the different possible vacua of the 
theory. 

These equations above are the same equations that can be derived from a one matrix 
model with potential V(X). This is, we have the matrix model 



/> 



[dX]exp(-Nfi- L tr (V(X))) (2.11) 

which is to be solved in the planar limit in the large N limit. We can write the loop 
operators 

< X k >= itr (X k ) (2.12) 
and the loop equations for the one matrix model 

< V'(X)X k >= fiJ2< X k - 1 - 1 >< X 1 > (2.13) 

i 

So it is immediate to see that they have structurally the same form as the ones given above 
(by setting Oj ~ fi < X 1 >). In the matrix model we have the normalization condition 
< X° >=< 1 >= 1, which identifies fi ~ < W a W a >= S, the gaugino condensate. It 
is the identification between loop equations of a matrix model and the chiral ring constraints 
that allow us to relate these two very different mathematical problems. 

Now, let us use the quantum moduli space technique to rederive the same loop equa- 
tions from a different point of view. The idea is to study the gauge theory with potential 



2.2, with one brane in the bulk, this is, with gauge group U(N + 1) x U(l). 

The classical moduli space is obtained by solving the F-term constraints. The con- 
straints for the A, B give us 

zAi - AiZ = B, t z - zBi = (2.14) 

and since z is a scalar, for A, B to get a non-zero vev and give us a brane in the bulk, we 
need that one of the eigenvalues of z is equal to z, and that A, B be eigenvectors of the 
matrix z with eigenvalue z. 

We can diagonalize z, and single out the (lxl) block matrices corresponding to the 
eigenvalue z. With this convention 



z = diag(z,zi,...z n ), B = {b u 0, . . . , 0) A { 



( Clj\ 



V / 



(2.15) 
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The other eigenvalues will be constrained by the equations 



Aifli - A 2 B 2 = V'(z) B X A X - B 2 A 2 = V'{z) (2.16) 
which result in V'(zi) = and 

Ml - b 2 a 2 = V'{z) (2.17) 

From here, the classical eigenvalues z% which were not singled out will be distributed only 
along the saddle points of the potential V. 

Now, to find the moduli space we need to write the constraints between the gauge 
invariant variables (with respect to the group U(N + 1)) given by z and the combinations 

a x b x + a 2 b 2 a x b x - a 2 b 2 
u = a x b 2 , v = a 2 b x ,w = , \l = (2.18) 

Clearly, the F-terms imply f2 = V'(z)/2. And we can also use the algebraic relation 

V'(z) 2 

w 2 = ft 2 + uv = — ^- + uv (2.19) 



to obtain the equations describing the moduli space given by the geometry 2.1 



The classical moduli space is given by the single equation in four variables w, u, v, z 

V'(z) 2 

w 2 -uv Y~ = (2.20) 

and it has branes located at the singularities of the geometry. Because there are singu- 
larities, it is possible to deform the moduli space and remove the singularities. These 
deformations are controlled by holomorphy. 

We should expect that the quantum deformed moduli space is of the form 



w 2 — uv = Quantum deformations = Q (2-21) 

where the right hand side is a polynomial involving only positive powers of the confining 
scale h-su(N+i)-> an d which is polynomial in the coefficients of P. This is the expected 
result from holomorphy [fTj], [48|, [l9| . 

The particular theory at hand has an SU{2) symmetry under which A\,A 2 and B±,B 2 
are doublets, while z, z are singlets, w, u, v together form a triplet, and O is a singlet. From 
these symmetries one can see that the quantum deformations will be independent of w, u, v, 
since terms of degree one in these variables that could appear in the polynomial Q above 
are not in a singlet representation of SU(2). Hence Q can only be given by a polynomial 
in z. The deformations should also be such that they are subleading with respect to the 
deformations induced by changing the couplings in the superpotential. If V'(z) is of degree 
n, and we change the functional form of V we can vary the coefficients of z n , z n+1 , . . . , z 2n 
independently of each other. The deformation should thus be of the form f(z), with / of 
degree smaller than n (]T2j |. 

It turns out that the number of parameters in the deformation / is equal to the number 
of roots of V . For each of these roots there is a conifold singularity, and each of these 



V'( 



,2 
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singularities can be deformed away by placing fractional branes at the given conifold, thus 



inducing a geometric transition which is of the same type as the one in [40]. We get in the 
end the same number of deformation parameters as there are singularities in the geometry. 

Now we want to see that these deformations can be related to the loop equations of 
the matrix model with potential V. 

For this, we need to know how to derive the quantum moduli space of the gauge field 
theory from the matrix model point of view. This has been previously argued for in 0]. 

We need to consider a multi-matrix model with the same field content as the gauge 
theory, and we take N — ► oo. However, we fix the matrix model so that the number of 
moduli stays fixed. This is, we will consider a matrix model with one probe brane in 
the bulk singled out. We do this with a matrix model potential given by the classical 
superpotential of the theory. 

{[dz][dz][dA)[dB])' exp(W) 

where the prime indicates that we leave some matrix fields unintegrated (exactly those that 
are massless perturbatively in the gauge theory). Since in the matrix model we do a large 
iV expansion, we are going to get to the results by analyzing the saddle point equations for 
the variables z, and what we called a^, 6j perturbatively above. Also, one of the eigenvalues 
of z in the gauge theory is singled out to be equal to z perturbatively. This eigenvalue will 
be called zq in the matrix model. At least naively, we can ignore the effect of the matrix 
probe on the large N matrix condensate because it would give rise to a 1/N effect. 

On going to a basis of eigenvalues for z, and including the Vandermonde determinant 
we find the following integral to perform in the large N limit: 

/n^o[^][^l 2 ][(^i2)dA 2 exp-{A>/i- 1 Ef = i(A J -^(^-^)+ (2-22) 
+AV 1 Ef=o V(\,) - N^ l V(~z) + Nh-^zq - 2)((oi)(6i) - (a 2 )(6 2 ))} (2.23) 

The above equation includes the Vandermode for all of the eigenvalues of z, not just the 
Vandermonde of the large N condensate. In the equations above A 2 is the Vadermonde 
determinant. The sum over j of A 3 Bj is over the color indices associated to the eigenval- 
ues Xi, . . . Xfj, but not to zq. This has been separated in the last line because we want 
01,61,02,62 to be moduli in the matrix model, so they cannot be massive at the saddle 
point. From this condition one of the eigenvalues of z is equal to z. 

In the equation above it is possible to integrate A 3 , Bj completely. This being a 
Gaussian integral over four coordinates of mass Xj — z gives us a measure term equal to 

^IWAi-*)- 2 

The logarithm of the Vandermonde determinant combined with 6 is then 

log(A 2 o) = 2 log(A 4 - Xj) - 2 £ log(A 4 - z) 

and it is the term in the effective action for the eigenvalues A and z. Now we want to solve 
for the saddle point of this setup in the large N limit. This will result in summing all of the 
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planar diagrams of massive fields in the above theory. Notice that the eigenvalues z and zq 
do not have an interaction between them, because we have not integrated out the massless 
modes of A,B. The saddle point equation for the zero components of A, B generic make 
z = zq in this situation. 

One sees that the saddle point equations for the eigenvalues Aj are the same as when 
we have the theory U(N) x U(0), because the contribution from zq cancels the contribution 
from z when they are equal. This seems accidental in the matrix model. The saddle point 
equation for the eigenvalues Aj is 

N»- V(A<) - 2 3— = (2.24) 



•3 



Now let w(X) = i X^i^o X^A" ^ e ^ ne resolvent of the matrix model. 
The saddle point equations for z = zq are identically equal to 

V'(z) - 2fiw(z) + 0161 - a 2 b 2 = (2.25) 

From here, we can use the same identity that led us to the classical moduli space to find 
that the quantum moduli space, as described by the matrix model, should be given by the 
equation 

w 2 -uv = -{V'(z) -2fiw(2)) 2 (2.26) 

Now, we compare the matrix model prediction to the prediction obtained by holomorphy 
arguments. It should be the case that two different functions of z are identical, these are 
as follows 

-{V'{z) - 2»w(z)f = \{V'{z)f + f(z) (2.27) 

where / is a polynomial. We obtain from these the following set of equations (after changing 
the normalization of / to / = —Afi^ 1 /) 

V'{z)w{z) - fiw(z) 2 = f(z) (2.28) 



which encode the full set of loop equations of the one matrix model [2.10 , when we realize 
that R(u) and w(z) are really the same object. 

To summarize: the quantum deformed moduli space of a supersymmetric gauge theory 
can encode loop equations for a matrix model. The above result for a one matrix model 
is not new and seems redundant, but the results in the following sections, using the same 
techniques are new. 

3. Superpotential deformations of J\f = 4 SUSY gauge theories 

The TV = 4 SUSY gauge theory is a supersymmetric field theory characterized by a gauge 
group G 2 and three fields X, Y, Z in the adjoint whose superpotential (up to normalization) 
is given by 

tr (XYZ - XZY) (3.1) 

2 The gauge group will be U(N) for the remainder of the paper. The results presented here can be 
generalized to other gauge groups which admit a large JV description 
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This theory is a super confer mal field theory which spontaneously breaks the conformal 
invariance at generic points in the moduli space, but if the vacuum is at the origin in moduli 
space (the gauge group is unbroken) then the theory has a vacuum with superconformal 
invariance. 

The gauge theory with N = 4 SUSY has a moduli space of conformal field theories 
parametrized by the gauge coupling of the theory g. If we consider deformations of the 
superpotential, this moduli space of TV = 1 superconformal field theories possesses a three 
parameter family of deformations which give rise to superconformal field theories [44]. 
These have (up to normalization) the following form of the superpotential 

W = tr {XYZ - qXZY + ^(Y 3 + Y 2 + Z 3 )) (3.2) 

and they also respect a Z3 symmetry which permutes X — > Y — » Z. When (3 = on top 
of the U(1)r symmetry, there is an additional non-anomalous U(l) 2 global symmetry of 
rotations of X, Y, Z by phases. 

It has been shown that some of these can have very interesting brane descriptions 
and AdS/CFT geometric duals [^, H, of the form of branes on C 3 /(Z n xZ„) or AdS*, x 
S 5 /(7* n x Z n ) for (3 = and q n = 1. These are geometric theories with branes on orbifold 
singularities, and there are three lines of singularities in codimension two that extend 
to infinity (one can think of them as the X, Y, Z axis). From here one finds a host of 
deformations of the geometry which serve to remove the codimension two singularities, 
leaving only codimension three singularities behind 3 . 

Since these examples are noncompact Calabi-Yau manifolds, one has to take care in 
defining the set of allowable deformations that one can study, as the total number of 
holomorphic functions describing the deformations is infinite. The deformations of the 
geometry modify the superpotential of the brane, so it is possible to analyze the geometry 
given a superpotential. 

For the rest of the paper we will consider only a special form of the general superpo- 
tential described above: 

W = tr(XYZ — qXZY + V(X)) (3.3) 

Where q n = 1 and 

V{x)= £ %i (3.4) 

i^O mod (n) 

is a polynomial potential. 

The first condition guarantees that there is a non-anomalous U(1)yz symmetry under 
which Y, Z have opposite charges, and is a remnant of the SU{4:) R-symmetry of the TV = 4 
gauge theory. This symmetry serves to simplify the problem. 



3 These codimension three singularities are usually of the conifold type, see the examples in [^, ^] 
4 One could also study by the methods presented here a situation where W = tr(XY Z — qXZY + V{X) + 

m(Z 2 + Y 2 ) + aZ + a'Y) but this example requires a lot more algebra to understand, the one cut solution 

has been studied in 
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The second point to notice is that when the Oi = we can use the other £7(1) charge 
combined with the R symmetry in such a way that the field X has different charge than 
Y, Z and we can also choose it in such a way that X, Y, Z have all positive charge as well 
as all of the coefficients in V(X), when we introduce them and assign quantum numbers 
to the a, so that we keep the global symmetry. Let us call this charge U{1)ri 

From this point of view, the deformation described above is a relevant deformation 
of the gauge field theory (when considering the grading of potential terms given by the 
charge), because if we consider the coefficients as ordinary c-numbers, the R' charge of the 
deformation is smaller than the charge of the original superpotential. 



The second condition given by equation 3.4 is more tricky to understand. From the 



string theory point of view it amounts to the deformations of the geometry being given 
exclusively by turning fields in the twisted sector of the orbifold which have charge zero 
under the U(1)yz symmetry. These guarantee that the dimension of the moduli space of a 
brane in the bulk remains the same as in the theory without the deformation, namely, the 
deformation can be understood geometrically and it does not produce an effective potential 
on the moduli space. In what follows we will argue what form the moduli space takes when 
we turn these deformations on. 

3.1 Calculating the classical moduli space of a probe brane. 

As a super symmetric problem, we need to solve the F-term constraints resulting from the 
classical superpotential as described above. These classical equations that determine the 
vacua are given by 

XY - qYX = 

ZX - qXZ = (3.5) 
YZ - qZY = -V'(X) = - a i xt 

i^O mod (n) 



and for a probe brane in the bulk the rank of the group TV turns out to be equal to n [26]. 

In particular if we think of X, Y, Z as formal symbols instead of matrices, we obtain 
the result that finding solutions to the equations above in terms of matrices is the same 
problem as finding representations of an associative algebra over the complex numbers, 
generated by X, Y, Z, and subject to the relations |3.5| |5|, q], which are exactly of dimension 
n. We are particularly interested in irreducible representations of the above algebra of 
dimension n, because this is the rank associated to a brane in the bulk. The irreducibility 
implies that the gauge group is broken to U(l) in the infrared, which is the gauge group 
associated to a single brane in the bulk. 

Finding these representations will let us calculate the classical moduli space of the 
theory and recover the Calabi-Yau geometry. The holomorphic functions on the Calabi- 
Yau geometry will be given by functions which are in the center of the above algebra 
H g, and the relations between these functions will be the relations that follow by 
virtue of imposing the equations |3.5| . In an irreducible representation all of these variables 
in the center can be shown to be proportional to the identity due to Schur's lemma. Hence 
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these are essentially gauge invariant, as they can be recovered completely from their trace. 
From this point of view, the relations in the center of the algebra are relations in the 
classical chiral ring of a one probe brane system. 

We need to solve the problem of finding the set of irreducible representations of di- 
mension n of the above algebra and the center of the algebra. To do this we need to solve 
for the F-terms above in dimension n. 

Solving the first equation is easy if we introduce the clock and shift matrices 
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(3.6) 



which satisfy QP = qPQ, Q n = P n = 1. These two matrices are linearly independent and 
generate the matrix ring of re x n matrices. A general re x n matrix M can be written in a 
unique form as a linear combination M = Y11j=o o-ijP l Q : ' ■ 

One can show that the first equation in [T^ is solved (generically) up to equivalence by 
X = xP and Y = yQ. This is the case when Y n = y n ^ and X n = x n ^ 0. To solve the 
second equation, the most general solution has to be of the form 

n-1 



Z = zQ- 1 P- 1 + Y,*iQ~ l P i ~ l 



(3.7) 



If we now substitute this expression in the third equation, we obtain that 

Y^VZiP^il - q l ) = -Y,^' 1 ^ 1 (3-8) 

Since the powers P l are linearly independent, we can solve for the z$ as follows 

Zi = -o^-V-Va - q*) (3.9) 

Thus we find that the generic irreducible representations of dimension re are parametrized 
by three variables x, y, z. However, one can show that various values of x, y, z can corre- 
spond to the same irreducible representation. The equivalent classes can be distinguished 



by evaluating the coordinates in the center of the algebra. Notice that the equation 3^ 
can only be solved if no term proportional to P^ 1 appears on the right hand side of the 
equation. This is exactly the restriction that was imposed on the form of the superpotential. 

We have found above a three-parameter family of irreducible representations. However, 
gauge invariance makes it possible to change x — > qx with some other modifications on y 
and z, so the variables above describe a cover of the moduli space. We need to evaluate 
gauge invariant combinations of x, y, z. 

It is easy to show that u = X n = x n is in the center of the algebra. From the 
representation above it is also easy to show that v = Y n = y n is proportional to the 
identity in the representation. Since Y, Z appear essentially on the same footing in the 
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equations 3.5 , then it follows that uu = Z n is also in the center; although it is difficult to 
write an expression for Z n in terms of the at, x,y,z. Finally, there is another variable in 
the center, which we will call t and whose form is given by t = XYZ + f(X), where / a 
polynomial. XYZ can be readily calculated in the representation above to give 



n-1 



XYZ 



xyz 



Z iV xpi 



xyz 



5>a 



1 x i P i = t - f{X) 



(3.10) 



i=i 



so that f(X) = - q l )~ l X l , and t = xyz. 

We see that in this particular matrix representation u, v, t are easy to calculate, and 
w is hard to calculate. Also, these are four variables determined algebraically by three 
parameters x, y, z, so they must satisfy some algebraic relation. 

In the particular case where V'(X) = 0, one easily finds that w = Z n = z n , and then 
the variables u,v,w,t satisfy the following equation 



uvw = t n 



(3.11) 



With the addition of V, the relation above is deformed. The importance of the calculation 
done above was to find the expression that described the variable t. 

The classically deformed relation above is relatively easy to calculate if we choose to 
write Z in a slightly different form, so that w can be read more easily 
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(3.12) 



\ Cn-l 0/ 

Then we find that Z n = J\Ch an d we can write equations for the Q in terms of t. Namely 

/Co \ 



XYZ = xy 



Ci 



V 



t-f(x) 



(3.13) 



CnJ 



which reduce to xyQ = t — f{q l x). Taking the product of all of these equations we find 

J[Ci = uvw = l[(t-f(q i x)) (3.14) 



x n y n 



Now, because the right hand side is invariant when we change x — > qx, the right hand 
side only depends on u = x n and t. So we find that the deformed geometry corresponding 
to the change in the superpotential generated by V is the following hypersurface in four 
complex variables u,v,w,t 

uvw = Q(t,u) (3.15) 
where Q is a polynomial in t and u of order n in t and of the same order as V 1 in u. 



-14- 



One can also notice that the polynomial Q(t, u) is not a generic polynomial in t, u of 
the given order, but it depends on fewer parameters (only the variables a,). 

When this happens, the geometry described above has some special properties that 
make it non-generic. In the above case, one can expect that the non-genericity of Q(t,u) 
is related to the existence of singularities for the curve in the u,t variables described by 
Q(t, u) = 0. The singularities of the CY geometry can be obtained by demanding that 
the partial derivative with respect to the variables u, v, w, t of the hypersurface equation 
vanish. 

The derivatives with respect to v, w give rise to the following equations 



uv 



uw 







(3.16) 



These can be solved for general u^Qhyv = w = Q. Under these conditions the CY 
has singularities at some values of t, u for which the curve Q(t, u) = is singular. These 
singularities are in general of the conifold type. These are easy to find: just look for 
repeated roots of t. These are values of x such that f(q r x) = f(x) for some r / 0. 

A second set of singularities is produced by setting u = 0. In this case it is easy to show 
that Q(t, 0) = t n , and for n/1 the derivative of Q with respect to t vanishes at t = 0. This 
produces a curve of singularities characterized by the equation vw = Q u (t, u)\ u= o. This 
curve degenerates to two of the three lines of singularities that are present in the orbifold 
geometry when Q u {t, 0) = 0. This is a curve of A n _i singularities. 

The singularities which will be important for us are the isolated singularities that we 
described first, and which are of the conifold type. At a conifold singularity there are 
usually two types of fractional branes present (even in the case where the singularities 
seem to have discrete torsion [||), which wrap the local cycle of zero geometric size with 
the two possible different orientations. 

The fractional branes in the above situation are obtained by taking a brane in the bulk 
to a singularity. The brane representation then becomes reducible and can be written as a 
direct sum of two sub-representations lim R = , R~ . Each of these sub-representations 
behaves in the classical low energy theory exactly as in the conifold field theory described 



in [41], so all of the adjoint fields are massive when we consider these representations. 

Let us write these representations for generic q. The idea is that the algebra of the 
F-terms will be a slightly modified SU q {2) algebra, so the solutions will look essentially 
like representations of SU{2). We can choose these to be of dimension k, where Y, Z are 
like ladder operators. 
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The commutation relations with X force X to have the following form 
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(3.17) 



X = xdiag(l, q, q 



- 1 ) 



(3.18) 
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Finally, the commutation relations of Y, Z give linear relations that overdetermine the 
UiCi- It is even possible to choose a gauge where yi = 1, so that only the Q are variables. 
Consistency of the equations gives a finite number of roots x that solve the problem, and 
the number of different roots is equal to the degree of V'(x). Irreducibility requires that 
all of the Q are different from zero. For q n = 1, if k > n one can show that the above 
representations are reducible, because the Q end up being periodic and a zero invariably 
shows up. 

The ones that survive for q n = 1 are exactly the ones that are associated to fractional 
branes in the geometry, and have rank less than n. If we evaluate t as defined in the bulk 
for the fractional brane representations, we find that t = f(q k ~ 1 x); and the condition on x 
to be a "root" are that f(x) = f(q k x). 

Notice that there is no gauge symmetry sending x — ► qx for the fractional branes. This 
means that a fractional brane contributes to traces of the form tr (X k ) for all k. These 
gauge invariant variables are not traces of elements of the center of the algebra. 

Also, at roots of unity for q, of the infinite number of representations of the slightly 
deformed SU q (2) we have for general q, only finitely many survive, and these are exactly 
the representations of dimension less than n. This is the property that will give a finite 
number of cuts later in the matrix model. 

The fractional branes can lead to isolated vacua of a U(N) theory for arbitrary N. 
These are exactly like the representations of SU(2) that are familiar from the J\f = 1* 
system. There are also fractional branes at the curve of singularities. These have a moduli 
space of vacua which is given by a connected n-fold cover of the curve. Given n branes in 
the curve of singularities one can move them so that they sit on top of each other in the 
CY geometry, but that each of them is in a different leaf of the cover. This configuration 
is produced by taking a brane in the bulk to the line of singularities. 

The description of the classical vacua of the theory for U(N) is given by splitting into 
irreducible representations of the algebra. 

R = ®NiRi 

and lead to a theory which in the infrared has a Yl U(Ni) gauge symmetry, with iV = 
Ni dim(Ri). If the Ri are branes in the bulk or in the curve of singularities, then one has 
a moduli space of vacua. Even if the R{ are fractional branes, R^, there can be a moduli 
space for these branes if they are accompanied by R™ p . 

The branes in the bulk behave in the low energy physics like Af = 4 field theory, 
so we can take them to a generic point in moduli space where they are away from the 
singularities, and also remove paired fractional branes to a generic point in moduli space in 
the bulk. After this is done, there are some fractional branes that can not be removed from 
the singularities by moving inside the classical moduli space. The number of branes stuck 
at each isolated singularity is constant in the moduli space (this is a classical statement). 

In this form we realize that the description of the moduli space of the U (N) theory 
comes in branches of different dimension, and the branches are characterized by how we 
choose to distribute fractional branes at the isolated singularities. 
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Among all of these possibilities, one can always find solutions where there are only 
fractional branes stuck at the isolated singularities. This is because there is at least one 
fractional brane with dim(i?Q,) = 1. In general we can have various isolated vacua where 
the fractional branes are distributed in some form at the singularities. These isolated vacua 
are exactly of the type that one believes one understands, because in the infrared one has 
pure gauge theory degrees of freedom. The SU(Ni) gauge groups should confine and the 
low energy theory reduce to ]^[ , just like in the J\f = 2 SYM softly broken to M = 1. 

Geometrically, if we place fractional branes at conifold singularities, these are expected 
to take the geometry through a geometric transition and deform the singularity away, if 
they have the same behavior as described in [[!(], j33| . This produces a new geometry where 
the quantum effects of the branes are taken into account. We will argue later, based on 
holomorphy that this should correspond to having a quantum corrected curve Q(t, u), with 
a more generic polynomial of degree n in t and s in u. 

We can now ask the question: how will we be able to see this geometry? The intuition 
is that if we use a probe brane, then the brane will have as it's moduli space the deformed 
geometry. This intuition is valid so long as we can ignore the back-reaction of the branes 
at the singularities to the probe. From the string theory point of view, this will be a large 
Ni limit. However, we will see that we can predict using matrix models that there is no 
back reaction even at finite Ni. 

The issue we should now consider, is that for every distribution of fractional branes we 
get a different deformation. In the string geometry this is because the geometric transition 
exchanges the algebraic cycle that the fractional branes are wrapping, to an S 3 with flux. 
The flux is given by the number of branes wrapping the cycle, and the flux tries to make 
the volume of the associated S 3 bigger in order to reduce the energy in the flux. Thus, we 
have to consider each vacuum configuration independently. 

In particular, the description of the quantum deformed moduli space of a probe brane 
has to be different for each branch of vacua. This is very similar to the discussion given 
in ||, where the order parameters that distinguish the vacua are equations which are 
only satisfied on some branches of the theory and not others. If we stop to think for a 
minute, the discussion we have done so far even at the classical level has this property. The 
equation of the CY geometry in the gauge theory assumes that we have split the classical 
representation of the vacua into irreducible representations, and the equation is only valid 
representation by representation. This is, we have made a choice of eigenvalues (we can 
think of block diagonal matrices as a generized form of eigenvalue), and the equations are 
valid eigenvalue by eigenvalue. 

However, the eigenvalues themselves are not gauge invariant, as we have to account for 
permutations of the eigenvalues. Only the symmetric functions of the eigenvalues are truly 
gauge invariant. However, the CY equations described above are not valid when we sum 
over eigenvalues. They are only valid eigenvalue by eigenvalue. To classically obtain the 
eigenvalues from the symmetric functions amounts to solving a polynomial and choosing 
a labeling for the roots. This extra label is the non-gauge invariant information that tells 
us that classically we are not studying universal equations on the chiral ring of the gauge 
theory, but rather equations in the chiral ring that are valid only on a particular branch of 
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the theory. 

This subtle point is very important. It means that the non-universal equations that 
determine the quantum moduli space does not only receive contributions from instantons, 
but it can be corrected by other strong dynamics effects, such as gaugino condensation. 

3.2 The quantum moduli space 

We have exhausted the discussion of the classical moduli space of the theory. Now we 
want to study the moduli space of one probe in the presence of fractional branes at the 
singularities, whereas in the previous section we discussed essentially the moduli space of 
one probe alone, plus the classical description of the branches of the gauge theory. We 
want to understand the quantum effects of the fractional branes on the geometry, and we 
will use a probe brane to test the shape of the geometry. 

Before we study these quantum deformations, one should make a field theory argument 
for the moduli space not being lifted, namely, for the absence of a potential on the moduli 
space. 

The argument one can make is semiclassical. If one examines the theory at infinity in 
moduli space, all the massive fields that connect the brane in the bulk to the fractional 
branes have a mass much higher than any confining scale. This means one should be able to 
integrate them out completely and reduce the problem to low energy effective field theory. 
The effective field theory on the probe brane is essentially TV = 4 SUSY gauge theory 
locally, and we do not expect that to change too much. The effective superpotential of the 
brane in the bulk should decay at infinity, because we are in a semiclassical regime. Now, 
the vev of the superpotential is controlled by holomorphy, and it vanishes at infinity. If 
the moduli space is lifted, then the superpotential must be a holomorphic function on the 
moduli space, and it should have poles or vanish exactly. The only possibility of having 
poles is if there are codimension one singularities on the classical moduli space. However, 
the moduli space described above is a Calabi-Yau geometry, which cannot have singularities 
in codimension one. This means that there is no effective potential generated on the moduli 
space, and the moduli space is not lifted. This does not prevent the moduli space from 
getting corrected. Indeed, our argument in the previous subsection suggests that it does, 
and when it does, the constraints that the deformed moduli space satisfies are controlled 
by holomorphy. This argument should be generic for moduli spaces of bulk-branes in a 
geometrically engineered theory, as in that case the moduli space of vacua of a brane in 
the bulk is always a CY geometry. 

To check the possible shape of the geometry we need to consider the U(1)r/ charge we 
defined. It is clear that when we keep the charges of the coefficients of V into account, the 
equation describing the moduli space is homogeneous for the U(1)ri charge. Since this is 
a non-anomalous symmetry of the gauge theory, this remains true in the quantum moduli 
space. A (partial) gaugino condensate <Sj will have the same charge as the variable t (this 
is the charge of the superpotential term). 

However, u can have a charge that is not commensurate with S, because we had some 
freedom in the definition of the U(1)ri charge. Also v,w are charged with respect to the 
U(1)yz symmetry, so they must appear together. 
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Holomorphy implies that all of the couplings, and S appear with positive powers in 
the deformed equation. Also, the moduli space is not expected to develop new branches, so 
it must also be polynomial in u,v,w,t, and it should reduce to the classical moduli space 
if we set the Si to zero. The R' charge of the S allow us to replace monomials where a t 
appears with an S; the coefficients depend on the branch of moduli space one is studying, 
but we can not change the coefficient of vw nor include higher powers of vw. This results 
in a deformed CY geometry where the equation that describes the geometry has the form 

uvw = Q(t,u) + 5(t,u) = Q(t,u) (3.19) 

where we have deformed the curve Q(t,u) to Q(t,u), making sure that the R' charge is 
conserved when we keep the information of the charge of the gaugino condensate. 

Now, let us take an example to illustrate the shape of the deformations, where we set 
n = 2 and s = 3. In this case V(x) = a\x + as/3x 3 , 2/(x) = a\x + a 3 x 3 = —2f(—x) = 
x(ai + a$u), and the curve is given by 

Q(t, u) = (t- x/2(a 1 + a 3 u))(t + x/2(ai + a 3 u)) = t 2 - u{a x + a 3 u) 2 /A (3.20) 

The singularity of Q(t, u) occurs at u = — a\/a 3 , t = 0. In principle the allowed deformation 
of the polynomial are given by 5(t,u) = r$tu + r\t + j3\u + (3q. The other coefficients, 
which are of higher degree, are determined uniquely by 0,1,0,3, the bare parameters of 
the superpotential. ro and t\ can be eliminated by a linear change in variables in t, so 
we drop them. In general we can always eliminate any term proportional to t n ~ l , so 
Q(t, u) = t n + t n ~ 2 1 (u) + t n - 3 2 (u) + .... 

Moreover, we have not placed branes at the curve of singularities located at u = 0, 
so we are not allowed to deform these singularities away. This forbids us from adding 
the term /3q. In the end the only allowed deformation is given just by (5\u. This is a one 
parameter deformation of the geometry. This is expected because we have only one conifold 
singularity to resolve. The counting of parameters that deform the geometry should be 
equal to the maximum number of conifold singularities that the geometry can have. When 
these conifolds collide they can generate more complicated singularities, but the counting 
of normalizable deformations should remain invariant. 

In the general case the polynomial Q(t,u) has no term added which does not contain 
a power of u, because we do not place fractional branes at the curve of ^4 n _i singularities, 
only at isolated singularities. This gives rise to constraints on the coefficients of the curve. 
The reader can verify that the deformed CY geometry still has a curve of singularities at 
t = u = with the above prescription of not including any deformation with just a different 
power of t which does not involve u as well. The general idea is that placing fractional 
branes at a given singularity only resolves that singularity, and leaves the other singularities 
intact. In fact this argument has been used in various papers |], |8| to study factorization 
properties of Seiberg-Witten curves. The geometrical idea is that one stabilizes the 3-cycles 
by flux, and in the absence of flux, the cycle shrinks to zero size. At least from this point 
of view one can characterize the configurations with no branes at some singularities: these 
are solutions where the curve Q(t, u) preserves these singularities, and is therefore singular. 
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We have presented a holomorphy argument to determine the deformed geometry as seen 
by a probe in the bulk, when we have placed fractional branes at the conifold singularities. 
The reader should be warned that the arguments given above assume that there is a way 
to define the variables u,v,w,t for a given configuration. Even classically, this depends 
on "removing" the information of the fractional branes to extract these variables. In the 
quantum theory this procedure gets corrections, so it becomes very tricky to define these 
in a sensible way. Even the partial gaugino condensates need a good definition in the field 
theory @. The underlying assumption of these ideas is that all of this can be done in a 
systematic way, and that there is a procedure to determine these variables in terms of vevs 
of elements of the chiral ring of the gauge theory. This has to be done in each particular 
branch of the theory separately, and with only one probe brane in the bulk. 

We will later show that as in the example of the previous section, the deformed geom- 
etry encodes relations for some generalized resolvents in a matrix model. This result can 
be used to find recursion relations which 'solve' the matrix model in the sense expressed 
in the introduction. 

3.3 The chiral ring of the gauge field theory 

We have argued what the shape of the deformed moduli space of vacua is, as it pertains 
to a single bulk D-brane. In some sense, this solves the structure of the chiral ring for a 
single brane in the bulk. 

As we have argued, the description of the classical moduli space depends on singling 
out eigenvalues and performing various operations with classical matrices. In the quantum 
theory all of these statements should be interpreted as operator equations in the chiral 
ring, so it is important to classify the chiral ring of the full gauge theory, not just of the 
"bulk" brane. 

The description is very similar to the one matrix model, except that now we have three 
matrices. Also, we should regard the rank of the matrices as an unknown, as we can have 
arbitrary fractional branes at the various singularities. 

The chiral ring is made of gauge invariant operators, and factorization is guaranteed 
by the cluster decomposition principle when we study chiral rings vevs in a supersymmetric 
vacuum. It is under these conditions that we actually believe that it is possible to relate 
the supersymmetric gauge theory to the associated matrix model. 

Following the setup of 0, the chiral ring is generated by traces of the following form 



where / is a polynomial in X, Y, Z. We have to remember that X, Y, Z are matrices, so 
the ordering inside / of the different letters X, Y, Z matters. The ordering of W a is not 
important, because adjoint action of W a on a set of letters is a total derivative. This allows 
us to choose an ordering where all the VF's are together. 



tr(f(X,Y,Z)) 
tr(f(X,Y,Z)W a ) 
tr(f(X,Y,Z)W a W a ) 



(3.21) 




(3.23) 
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Also, no more than two W a per trace are allowed because these terms are a total chiral 
D derivative 0. 

The chiral ring generators as presented above are independent of N, and of the specific 
details in the superpotential. The generalized Konishi anomaly equations for variations of 
the fields will depend on the specific form of the superpotential, and these are the ones 
that can be related to a matrix model. 

Classically, the relations in the chiral ring are given by solving the F-term constraints 
of the superpotential. These can be written in a gauge invariant way in the following form 

dW 

tv(—f(X,Y,Z)) = (3.24) 

for / any word in the fields, and <j> any of X, Y, Z. These are essentially the same equations 
as W = 0, since it implies that no matter what matrix we build out of X, Y, Z, the trace 
is zero. 

If we include quantum corrections, the equations of motion take the form 

tr (W f) = Quantum corrections (3.25) 

The quantum corrections involve the (partial) gaugino condensate (s) of the (fractional 
branes) gauge theory, and we should interpret the above equation as a part of a generalized 
Konishi anomaly, where / can include the W a W a . 

A few points are worth noticing: when we fix the superpotential as in the previous 
subsection, branes in the bulk do not contribute classically to one point functions of the 
form 

tr(X k ),tr(Y k ),tr(Z k ) (3.26) 

for k 7^ mod (n), nor to any word containing a gaugino field 5 . This can be seen because 
the classical solutions of branes in the bulk give contributions to the above as sums of the 
form 

£<f fe = 

i 

In the orbifold theory, all of these words correspond to twisted fields. The classical 
string theory statement is that a brane in the bulk does not carry charge in the twisted 
sector of the orbifold, so it can not generate vevs for twisted fields. 

Fractional branes, on the other hand, do couple to the twisted sector. It follows that 
they generate vevs for twisted sector words. We should read this statement in the following 
way: classically, there is information that depends only on the fractional brane structure. 

Moreover, we are interested in isolating the fractional branes that are stuck to the 
codimension three singularity, in a situation where there is no moduli space of vacua. In 
the gauge theory the low energy physics reduces to pure gauge theory Y\U(Ni), and we 
can argue that the SU(Ni) confine and give rise to some partial gaugino condensate. 

a Classically the gaugino field has to be zero for Lorentz invariance, and this implies that the gaugino 
condensate classically is also zero 
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Guaranteeing that there are no branes in the codimension two singularities of the 
geometry is a simple matter. These contribute to the vevs tr (Y k ) and tr (Z k ) for all k. 
Also, branes in the bulk can contribute to tr (Y k ),ti (Z k ),tv (X k ) for k multiples of n. 

This means that classically these can be considered variables that we can set to any 
value we want. Notice that tr (X n ) and tr (X 2n ) are uncorrelated if we have more than 
one brane in the bulk. In particular, we can always consider enough branes in the bulk so 
that these variables are allowed to have any value. 

The same is not true for the other powers of tr (X k ). There are only a finite number 
of values of X that are allowed for a fractional brane, so these variables are related to 
each other in the classical theory; there will be recursion relations that determine these 
variables. The data necessary to seed this information to the classical recursion relations 
is exactly how many fractional branes there are at each singularity, and which cycle they 
are wrapping. 

One can wonder if in the quantum theory the vevs of tr (X k ) and tr (X k W Q W a ) for 
k not a multiple of n can get contributions from branes in the bulk. We will argue that 
for the second ones, which have a simple matrix model interpretation this is not so. A 
hand waving argument would say that the gauge groups that have gaugino condensation 
are just the fractional branes and not the branes in the bulk, so they contribute to these 
one point functions. In principle, if these were classical matrices, one could make the same 
statement about any trace containing W a W a . However, the quantum vacuum is not given 
by classical matrices, and one should consider only twisted operators. 

At infinity the branes in the bulk should be semiclassical and decouple, so their ef- 
fective gaugino condensate should vanish. Thus the contribution to tr (X k W 2 ) should be 
suppressed by a large mass scale. At infinity, their contribution is zero, and tr (X k W 2 ), 
being holomorphic on the moduli space, is constant by the same argument that made the 
effective superpotential on the moduli space equal to zero. This suggests that these one 
point functions are determined by the fractional branes alone. 

Since these are related to tr (X k ), and classically the twisted tr (X k ) satisfy recursion 
relations (they are determined alone by the fractional branes), it suggests that the same is 
true for tr (X k W 2 ). The description above is suggestive. It is not meant to be a formal 
argument. We will use matrix models to make this intuition precise. 

3.4 Relations to a matrix model 

We have found that we can write the most general form of the deformation of the Calabi- 
Yau geometry which is compatible with holomorphy. This information in the gauge theory 
encodes non-perturbative dynamics that results from gaugino condensation on the gauge 
group of the fractional branes. All of these effects should be computable from a matrix 
model. What we have described is the moduli space of a probe brane. The matrix model 
calculation proceeds in the same way we worked in section We isolate the block of 
matrices associated to the probe brane. This is done by (perturbatively) identifying a set 
of eigenvalues for X ( an n x n block matrix) to correspond to a bulk brane, where the 
eigenvalues are related to each other in such a way that they produce massless field in the 
classical field theory. 
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In the end, we have to consider the three matrix model whose potential is given by 

W(X, Y, Z) = N^r 1 [tr (XYZ - qXZY) + tr V(X)] (3.27) 

with n correlated eigenvalues which correspond to a single brane in the bulk, and take the 
large N limit. 

We are after the effective potential for the set of block matrices in the bulk. Notice 
that we can choose the gauge where X is diagonal, and we can isolate the block of matrices 
which correspond to the eigenvalues in the bulk. 

The idea now is to integrate the off-diagonal elements that connect the eigenvalues of 
the bulk to the eigenvalues of the fractional branes. The procedure is to integrate these 
fields in the matrix model, where we take a large iV limit for the matrices and we will 
obtain an effective potential for Xb u ik,Ybuik, Zbulk with the effects of the fractional branes 
included. 

This should be seen in the same philosophy of the proof of the matrix model conjec- 
ture as given in |L5|]. The new ingredient is that perturbatively massless fields can not be 
integrated out. We also have to remember that since the classical moduli space has singu- 
larities, at each singularity new massless fields appear. The nxn block of matrices contains 
all of these fields, so we need to keep them unintegrated to understand what happened to 
the singularities. 

To compute the quantum corrected moduli space then we need to vary the effective 
potential for Xf, u ik^buik, ^buik as n x n matrices and set the variations equal to zero. We 
then need to evaluate what are the quantum corrections to the form of u,v,w,t for the 
bulk matrices and find the new relations between them. This is technically a little bit more 
elaborate than the situation in section |2[ It is useful to keep in mind that the brane in the 
bulk is auxiliary: it is giving us information about the matrix model of the branes stuck 
at the singularities. We want to know exactly what information it is encoding. 

To do the first part, we divide the matrices as follows 

fX b \ 



V 

where we have taken care to write the matrix X in a block diagonal form, and where the 
eigenvalues attached to the singularities are singled out. We keep the matrix of the bulk 
Xb in general form. In the matrix model, with Xb unintegrated, requiring that the Xj,f = 
is a gauge choice. This introduces ghosts [JI6[j, whose measure will give the Vandermonde 
determinant. The subindices are labeled as a mnemonic: b for bulk, / for fractional. 
Integrating over the orbit of the gauge group introduces a factor of the volume of the orbit, 
which is proportional to the Vandermode determinant. This includes both the eigenvalues 
at the singularities and the block matrix Xb, even if we choose not to diagonalize Xb, there 
is still a contribution from how the eigenvalues A 'repel' the block Xb- 





(3.28) 
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When we take the logarithm of this measure term, it gives us a contribution to the 
effective potential which is given by 



2 ^ tr 6 log(X b - Ai) + - Xj) 

i i^j 



(3.29) 



We have only written the terms involving X b , as this is the information that relates to the 
brane in the bulk 6 . 

We now set X b f = everywhere in W(X,Y, Z) and remember to keep track of the 
Vandermonde determinant. 

It is easy to see that after this substitution Y b f,Yf b and Zf b , Z b f appear quadratically 
in the equation, and can be integrated completely, separating the matrices from the bulk 
from the matrices associated to the fractional branes, except for quantum effects. These 
two off-diagonal integrations contribute another effective term to the potential of X b alone, 
which is given by 



J> b {\og{X b - qh) - ^log(X 6 - q- x \)) 



(3.30) 



Again, we obtain a formal power series in X b which encodes loop operators of the fractional 
brane matrices. 

The effective potential for X b , Y b , Z b in the matrix model is then given by 



W = N/1- 1 (tr (X b Y b Z b - qX b Z b Y b ) + V(X b )) 



E 



2log(X b - Aj) - log(X b - qXi) - ^ ^og(X b - q^X,) 



(3.31) 
(3.32) 



Now, we have to assume that the distribution of the Aj is known, so that we can think of 
the terms involving logs as a formal function of X b . Notice that this formal function has 
a power series expansion in X b which has the same characteristics as V(X b ): no powers 
of X b which are multiples of n appear. The effect is the same as if we would have done a 
change in the potential 



N ^ 



2log{X - Xi) - log(X - qXi 



E 

i 



log(X - q^Xi 



so the moduli space has essentially the same form as we discussed in the previous section, 
except that the new Q now has information about the distribution of eigenvalues Aj. Notice 
that the formal form of the solution gives us a curve of the following form 

UVW = Q matrix (t,u) 

where in principle Q ma trix(t,u) is a formal power series in u~ l . 



"The saddle point equations for the eigenvalues A are discussed in the appendix. 
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The new calculation of V(X)' gives the following result 



V'(X) = V'(X)+ £ 



T2— - - 

Y X-Xi X-qX t X-q-^Xi 



(3.33) 



n-1 



X n ~i- l X j ; 



where u = x n = X n is still in the center of the modified algebra. 

We can define the following n — 1 generating functions of moments of the eigenvalues 

Xi, 

1 X J °° s \nk+j ^ 

< 3 - 35 > 

i ' 1 k=0 ' 

and write the (matrix version of the) quantum moduli space in terms of these generating 
functions. 

The quantum version of f(x) is 

n-1 

f(X) = f(X) + £(1 - ^)X^i?» (3.36) 
j=i 

Now we can compare both descriptions of the quantum moduli space: the one we 
obtained using holomorphy arguments, and the one where we assumed that we had the 
knowledge of the Rj(u). The main assumption in this paper is that both constructions 
are giving exactly the same moduli space. Indeed, the one difference between Q(t, u) and 
Qmatrix{t, u) is that the first one is a polynomial, and the second one is only a formal power 
series. 

The equality of the moduli spaces means that we can write equations that the Rj (u) 
must satisfy: one for each power of t, so that the matrix model gives also a polynomial. 
There are n — 1 unknown functions Rj(u), and there are n — 1 non-trivial relations, one 
for each power of t starting at t n ~ 2 . The equation for t n ~ l is satisfied trivially because the 
product defining Q is over all the values of f(q l X), and / does not contain integral powers 
of X n . This means the solution is algebraic. 

The other equations are given by expanding both forms of Q(t,u) and comparing 
coefficients. We will call these the Quantum relations of the matrix model. The first two 
are explicitly written bellow: 



1 (u)= y £f(q i x)ftfx) (3.37) 

i<j 

-o 2 (u) = 2 ftq i x)f(q i x)f(q k x) (3-38) 

i<j<k 

The right hand side is given by explicit power series of in terms of the one point 
functions < X k >, and all of the coefficients except the first few (the ones that have 
positive powers of u) vanish. This gives us constraints that the one point functions must 
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satisfy, and one can notice that they give rise to recursion relations for the moments of the 
eigenvalues < X k >, so that all but a finite number of them are determined in terms of 
initial conditions. The polynomials Oi serve to encode these initial conditions. 
In more detail, we can decompose / in terms of powers of x modulo n, 

f j (x) = (l-qi)- 1 Yl aixi ( 3 - 39 ) 

i=j mod (n) 

Oi(u) = Y^jUl l .v)p:H , .v) (3.40) 

i<j 

= J2l[(f j (l-q j r 1 + Kl-Q- j )Rn- j (u)xi) (3.41) 
j 

(fn-j(l ~ q-*)- 1 + (1 - q j )fiR j (u)x n ^) (3.42) 
(q> +q 2j + --- + g ( n - 1 ^')(n-l)] (3.43) 

The sum appearing in the above equation can be evaluated using l+qi +q 2 i + ■ ■ ■ -\-q^ n ~ 1 )^ = 
0. We get 

1 (u) + (n - 1)^(2 - <j - q~>)-\f 3 {x)f n ^{x)) = (3.44) 



li(n - 1) 



(3.45) 



3 

Now, the left hand side of the above equation is a polynomial in u, while the right hand 
side is a power series in u^ 1 . Equality means that the coefficients of u~ k of the right hand 
side must vanish identically. 

After a little bit of algebra collecting the powers of u^ 1 this gives us the following 
relations 

Y, a i< A >= — (3-46) 

i 

Now, we will show that these equations can be derived from manipulating loop equa- 
tions in the associated matrix model. We can ask why these relations look so different from 
each other for different values of n, and what is special about q being a root of unity in the 
matrix model. This is what we will analyze in the next section. Notice that the concept 
of holomorphy has allowed us to propose some polynomial relations for twisted one point 
functions in the matrix model, but at this moment this should be interpreted as a guess, 
since there is no proof that the construction of the moduli space based on matrix models 
necessarily gives the right answer. 

4. Obtaining the Quantum relations from loop equations for a matrix 
model 

Consider the three matrix model whose potential is given by 

W{X, Y, Z) = iV/i _1 [tr (XYZ - qXZY) + tr V{X)} (4.1) 
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with 

V(x)= £ (4.2) 

i^O mod (n) 

a polynomial, and q n = 1 a primitive n-th root of unity. 

Because the matrix model is given formally by the superpotential of a supersymmetric 
field theory, we will classify the one point functions < f(4>) > by their quantum numbers 
under the Z n x Z n quantum symmetry of the corresponding chiral operator in the quantum 
field theory. A word which has charge zero under the quantum symmetry will be called 
untwisted, and a word which transforms non-trivially under the quantum symmetry will be 
called twisted. This nomenclature is borrowed from the fact that the theory with V(x) = 
results from the field theory at a special orbifold singularity. We will also assign Y and Z 
charges ±1, so that we will also use this U{\)yz charge to classify words. 

Since this second charge is conserved by the potential and the measure, we can consis- 
tently set to zero the one point function of any word whose charge is different from zero. 
This constraint at the level of the classical field theory implies that there are no fractional 
branes on the singular curves of the CY, and is a simplifying assumption to reduce the 
amount of algebra. We want to impose this constraint as well on the matrix model, since 
we are interested in situations where the matrix model has an isolated vacuum about which 
all of the modes are massive. 

We will use throughout the following simplifying notation convention for the one point 
functions of words in the matrix model 

< m >=< ^tr (f {</>)) > (4.3) 



Now we will proceed to derive the equations 3.37 from the loop equations associated 
to this matrix model. 

The universal loop equations can be written by considering an infinitesimal change 
of variables 5(f) 1 = f(4>) (one assumes this is a monomial written as some ordering of the 
matrix variables, i. e. a word in the matrices), and realizing that the following expression 
is the integral of a total derivative and therefore vanishes 



[dfld^SP)* eM-W) = J [#]tr {d+iSP) e W (-W) = (4.4) 
The loop equations will have the universal form 

Ea</i><A>=<^> (4-5) 

where the sum is over all fx, fi monomials such that f\<f) % fa = this is, from splitting 

the word / into two words connected by the letter (jj 1 . The above result is using factorization 
of the large N limit of the matrix model. 

Because we set all words with charge different than zero to zero, we need to look only 
at variations of the matrices X, Y, Z which have the same charge as the object we are 
varying. Thus we will look for 5X = f(X, Y, Z), where / has the same number of Z and Y 
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matrices. Similarly 5Y = g(X, Y, Z) where the word g has one more field of type Y than 
of type Z. 

There is also the U(1)ri to consider. It's existence gives a grading to the set of words. 
One can systematically explore the set of possible monomial variations by exploiting the 
U(1)ri symmetry and take variations with increasing U(1)ri charge. 

4.1 The first equation 

We will study the loop equations inspired in the solution of the two matrix problem given 



in [51|. The main idea is to write the loop equations in such a way that one can learn 
information on complicated words in terms of simpler words. 

Consider the following matrix variations, of minimal number of Y, Z pairs 

6Y = x m YX k ~ m (4.6) 

The loop equations associated to this variation take the simple form 

H < X m >< X k - m >=< X m+1 YX k - m Z - qX m YX k - m+1 Z > (4.7) 

On the right hand side there are k + 2 independent possible orderings of the fields, while 
we only have fe + 1 different equations. These one point functions on the right have exactly 
one pair of Y, Z matrices. Variations of Z of this type will produce the same equations as 
above, so they do not lead to new loop equations. 

Also, one can consider the variation of X given by SX = X k+1 , from which we get 

k 

H^2<X m >< X k - m > - < V'{X)X k+1 >=< X k+1 YZ > -q < X k+1 ZY > (4.8) 

m=0 



which gives an extra equation for some of the same one point functions in equation [4.7 . 

In total we have k + 2 equations for the k + 2 variables < x m YX k+1 ~ m Z >, so in 
a generic case we could say that the one point functions of < x m YX k+1 ~ m Z > are all 
determined from the one point functions of < X m >. Considering the f/(l)#/ grading, the 
words involving Y, Z have higher degree than those with only X. 

Notice that one can use the equations to shift Z to the right in the word < x m Y X k+1 ~ m Z >, 
and then eventually one comes back to the original configuration with a phase, so in this 
case one can solve for any of the one point functions in the right hand side, as long as the 
phase is not zero. 

However, there are situations in which these equations are not linearly independent. 
This will happen when k is a multiple of n, and then on transporting Z around the 
loop the total accumulated phase is q k = 1. From here the equations for the vevs 
< x m Y X k+1 ~ m Z > are linearly dependent. This means that the one point function 
jmj/jfc+i-m^ > i s not determined in this case. Notice that this is exactly the case 
when the word is untwisted 7 . 



7 From the point of view of the classical SUSY gauge field theory, this function in the chiral ring receives 
contributions from the moduli space of a probe brane, so it can take any value. 



- 28 - 



Consistency of the loop equations then imposes a constraint between the elements of 
the left hand side. This constraint is given by summing 

J^/V < X m >< X k - m > = ^ q - m < X m+1 YX k - m Z - qX m YX k - m+1 Z > 

= q- k < X k+1 YZ > -q < YX k+1 Z > (4.9) 

= < X k+1 YZ > -q < YX k+1 Z > (4.10) 

k 

= n < x m >< x k -' m >-< V\X)X k+1 > 
and this equation can be rewritten after some manipulations as follows 

< V'(X)X k+1 >=\Y1 ^ 2 -9 m - < X ™ >< xk ~ m > ( 4 - n ) 



which is exactly the first of the relations that appeared in |3.46 



As the reader can verify, the equation above is written only in terms of twisted one- 
point functions, as {2—q m — q~ m ) = whenever m is a multiple of n. Since V is a polynomial 
(lets say of degree s), the above equation can be seen as a recursion relation that determines 
the highest degree component of V'(X)X k+1 , namely the one point function of < X k+S > 
(in the case where k = rn), in terms of one point functions of lower degree. The first such 
component that can be determined is the one corresponding to < X s >, and then all other 
components modulo n. 

Now, we will solve the constraints for the one point functions < X k+l YZ > for k not 
a multiple of n, because we will need them for future calculations. This is easily seen to be 



given by solving 4.6 when k is not a multiple of n. In this case there is no equality between 
the second and third lines, and instead we get two linearly independent equations with two 
unknowns. We get that 

(q- k - 1) < X k+1 YZ > = ^fiq m <X m X X k ~ m > 

k 

-fi, < xm >< xk ~ m > + < V'(X)X k+1 > (4.12) 

m=0 

= fiY(i m - x ) < xm >< xk ~ m > + < y\ x ) xk+l ^4.13) 

Notice that on the right hand side of the third line we do get untwisted one point functions, 
so these can in principle receive contributions from branes in the bulk. This happens when 
k — m is a multiple of n, or from V'(X)X k+l when certain terms in the potential are 
non-zero. 

4.2 The second Quantum equation 

The next level of difficulty is obtained by taking words of next to minimal numbers of 
pairs of Y, Z fields. For example 5Y = YX m ZX n Y, and all possible orderings of such 
monomials. Although it looks much harder to solve this linear system, we do not need to 
study the most general variations of Y of this type. The main effect of the loop equations 
is to tell us what happens when we permute two of the symbols in a word in terms of 
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one point functions of lower numbers of pairs of Y, Z. Thus any one point function can 
be determined from any other permutation by a series of steps. We only need to worry 
about the possibility of the system of equations not being linearly independent, just as we 
saw previously. This can happen when the words we study are untwisted, and then one 
produces a constraint between words of smaller degree. 

To check this redundancy, we just need to take an initial word < X k+2 YY Z Z > 
and cycle the Z around the X and Y. We will return with the same phase to the initial 
configurations when A: is a multiple of n. This is because passing a Z through an X gives a 
factor of q, and passing a Z past a Y gives a factor of q^ 1 , so the total phase from cycling 
Z around the trace is q k+2 q~ 2 = q k . 

We only need to look at the variations SY = X m Y 2 Z X k+1 ~ m , and their corresponding 
loop equations 

fi < X m X YZX k - m > = < X m+1 Y 2 ZX k+1 - m Z > -q < X m Y 2 ZX k+2 - m Z > (4.14) 

where we have set the charged one point functions to zero. 

We also need two equations to cycle the Z past the Y. These are provided by the 
variations SX = YZX k+2 and SX = ZX k+2 Y . These two loop equations will look as 
follows 8 

< YZYZX k+2 > -q < ZYYZX k+2 > = 
nY,< Yzxm >< X k+1 - m > - < YZX k+2 V'{X) > 

rn 

< YYZZX k+2 > -q < YZYZX k+2 > = 

- < ZX k+2 YV'(X) > (4.15) 



Summing 4.14 we get, in the case where A; is a multiple of n that 
< x k+2 Y 2 ZX k+l-m Z > _ q 2 < Y 2 zx k+2 Z > = 

V Emio q k+l - m < X m >< YZX k - m > (4.16) 



which is also the sum of the equations in 4.15 



< x k+2 Y 2 Zx k+l-m Z > _ q 2 < Y ^ Z x k+2 Z >= 

^E m < YZX k+1 - m ><X m > 
-q < YZX k+2 V'{X) > - < ZX k+2 YV'(X) > (4.17) 

We can rewrite the constraint as follows 

p^ m (l - q- m ) < X m >< YZX k+l - m >= 

q < YZX k+2 V'(X) > + < ZX k+2 YV'(X) > (4.18) 

Again, we notice that the one point functions of < X m > for m a multiple of n do not 
contribute, because q m = 1. This means the terms proportional to /i only depends on 



3 We write only the equations with < X k Z >= for all k 
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twisted one point functions. Notice that the right hand side also only depends on twisted 
one point functions, at least as far as the powers of X that are excluded from V{X). Notice 
that the constraint above is not in the form which we want, which depends only on one 
point functions of the type < X m >. However, we can use the results from the previous 
subsection to express all of the one point functions appearing in the above equations in 
terms of one point functions of < X m >, possibly including the ones where m is a multiple 
of n, which we have no way of determining from our previous arguments. The equations 
obtained in this manner will contain terms of order /i 2 (triple trace), of order fj, (double 
trace), and of order fi° (single trace). We still have to check that after we do all of these 
substitutions, the one point functions appearing in the above equation are all twisted, and 
that the equation coincides with what we have shown already. This would lend support to 
the existence of a recursion relation determining only the twisted powers of X. 

At this point we need to substitute explicit forms for V' (X) in terms of it's polynomial 
expansion, because the way we will substitute in equation |4.12 depends on what particular 
powers of X are available. We also need to do an extra step in equation 4.18 and move the 
Z past some of the X in the one point function of < ZX k+2 YV'(X) >. 



Consider expanding V'(X) = ^ 



^0 mod (n) 



a>iX l . It is clear that we will get one 



point functions of monomials of the type < YX l ~ l ZX k+2 >, in the case where A; is a 
multiple of n. The following equality is easy to show for these one point functions by 
summing terms of the equations 4.7. 



< YZX i+k+1 > 



y- 1 < YX^ZX^ 2 >- 



i-2 
m=0 



This can be used to rewrite 



< YX l - l ZX k+2 >= q l ~ 



< YZX i+k+1 > -n <T < X m >< X 



fc+i+l— m 



> 



m=0 



(4.19) 



(4.20) 



From here we can check that the coefficient of ai < YZX k+1+l > in equations 4.15 coming 
from the potential is given by q + q % ~ 1 . 

We will now only check that the terms of order which are untwisted and propor- 
tional to the second power of the coefficients in the potential (a) 2 are canceled. For this, 
we need to use equation 4.12 and only keep terms that involve two powers of the a. 



l V'{X) > +0(n) 



Indeed, we get the sum 

zZil + 2 1_ >i < YZX k+i+1 >) = ^(q + q 1 -*)^-^ - l)" 1 ^ < X k+i+1 ^ 

i i 

= z2(i + - iH%- < * k+l+] > +o{n) (4.2i) 

Now we want to consider the case where i + j is a multiple of n, and we are summing two 
different terms with the coefficient OjOj < X k+t+ i >. Also, remember that k is a multiple 
of n, so we can drop it from the powers of q. The coefficient for the sum of these two terms 
is equal to 

q{l + - I)" 1 + g(l + q-i)(q-i - l)" 1 (4.22) 
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where we substitute i = —j mod n, so that the numerator of the expression is equal to 

q(l + q^q" 1 - 1) + q(l + g"*)^ - 1) = q{q~ i - q l + q i - q^) = (4.23) 

This shows that there are no untwisted powers of X appearing in the one point func- 
tions of this second constraint which are of second order in the aj. As these are also the 
terms that contain the highest power of X that is available. This shows that the above 
equations will determine higher twisted one-point functions of X. 

We still have not shown that in the double trace and triple trace operators one can 
eliminate all the untwisted powers of X. Since we have already found relations using 
the matrix theory calculation of the moduli space, and moreover we have shown that the 
constraints give the right result for the first quantum equation plus a consistency condition 
for the second equation, we will leave this lengthy algebraic calculation as an exercise to 
the reader. 

The methodology for determining the recursion relations should now be evident to the 
reader. One has to add more and more pairs of Y, Z operators, and consider variations of 
the form 

SY ~ X k+s - m Y s+1 Z s X m (4.24) 

which allow us to circulate one Z past all the X. We also need to consider the variations 
of X that allow us to flip Z past the Y in words with charge 0. Again, when the words are 
untwisted, the words will be undetermined classically and the loop equations are linearly 
dependent. Consistency imposes constraints that will involve twisted words with smaller 
numbers of pairs of Y, Z. Since these words are twisted, they will belong to the set of 
relations that can be solved in terms of smaller numbers of pairs of Y, Z variables. One 
should be able to show that after performing some amount of algebra one can reduce the 
system so that there are only one point functions for powers of X which are untwisted. 
The calculations at each succeeding step become more involved 9 , but this is a recursive 
algorithm to get the constraints. 

Summarizing, we see that the presence of a moduli space of vacua has given us a 
degenerate system of linear equations for the one point functions < X K Y M Z M ~ 1 X K> Z >. 
For generic values of q the associated linear equations determine these one point functions 
completely in terms of functions of lower degree. However, one does not get constraints that 
involve one point functions of X alone. These should be considered as unknown parameters 
determined by the "infinite number of cuts" of the matrix model. For the special values 
of q that we studied, this does not happen, and the linear system is degenerate. For this 
linear system to have a solution, the terms of lower degree must satisfy a constraint. It 
is exactly these constraints that encode the shape of quantum moduli space of vacua of 
the associated SUSY gauge theory. Moreover, the presence of only a finite number of 
singularities in moduli space guarantee that for certain generalized resolvents, one only has 
a finite number of cuts. Also, due to the special form of the potential, in the particular 

9 It would be interesting to find a more effective method to calculate these relations in the matrix model 
to give a more compelling proof, say in terms of eigenvalue distributions for X. The author has not found 
a better way to understand it however. 
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case we are studying, the solution effectively reduces to a one matrix model whose solution 
is determined by a Curve Qquatumiti^)- 

5. Partial gaugino condensates 

We have now described a partial solution of the three matrix problem, and we have found 
that some interesting generalized resolvents (the Rj(u)) can be obtained from the curve 
Q(t,u). 

In particular, it is interesting to ask how many eigenvalues in the matrix model are 
located around each classical saddle point of the theory. In the gauge theory this will be 
the calculation that determines the gaugino condensates from the geometry. 

Now we want to relate the Rj (u) to the problem of the distribution of eigenvalues of 
X. In particular, Rj(u) only has information about the distribution of fractional branes, so 
it only contains information about the eigenvalues of X associated to the fractional branes. 

Let us look at the definition of R n -\{u). It is given by 



i 1 k=0 



N u - A™ u 



Classically, this function has poles at the eigenvalues A™. It is more convenient to introduce 
a cover of the u plane where u = 7™, and then write 

R; ( 7 ") = iy — (5.2) 

i 1 

In particular, if we choose j = n — 1, the residue of the poles associated to a particular 
eigenvalue 7 = Aj is 1/nN. 

However, we have to take into account that classical saddle points for the eigenvalues 
associated to fractional branes R+ come located also at qx, q 2 x, . . . q r ~ 1 x where r is the 
dimension of the associated representation. And there is a conjugate representation R~a 
that appears as well for the other set of eigenvalues q r x, . . . q n ~ 1 x. 

So if we want to count how many eigenvalues N a in the matrix model are located 
around a classical saddle point which has an eigenvalue located at xq, we want to do the 
following contour integral 

- £ ^(7*^7 = - v^d + + ■ ■ ■ <r +1 > (5.3) 

where C is a contour over the eigenvalue plane which contains Xq. Notice that the above 
expression only depends on N+ — N~ which is the total number of fractional branes that 
can not be removed away from the singularity classically. If one studies the Klebanov 
Strassler system, this would count exactly the difference between the ranks of the gauge 
groups, which is an invariant under the duality cascade. 

Giving specific values to N a /N should determine the curve in the matrix model, by 
asking which curve has the right residues for a given value of \i. 
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In the gauge theory, the generalized resolvent R n -\ in the matrix model is related to 
the generalized resolvent for the chiral ring given by 

XiW 2 

K^ n )=tv(^-^) (5.4) 

The gauge invariant definition of the partial gaugino condensate Si associated to the 
fractional branes at some conifold singularity should be given by the following expression 

S a = n ; r tt I 7e n _i(7 n )d7 (5.5) 

2m(l + q- 1 + . . . q- r+1 U ' 



c 



where the contour is around one of the eigenvalues of X at the given singularity, and S a is 
the partial gaugino condensate of the full collection of fractional branes at the singularity, 
once we get to the bottom of the duality cascade of the local conifold. This expression 
should be compared with the recent setup of Dijkgraaf and Vafa based on supergroups 



1 2C ] , where everything only depends on the "difference of the gaugino condensates" . That 
statement is equivalent to branes in the bulk not contributing to the deformations of the 
geometry, as seen also in [||. 

From here, again, it should follow that one should be able to determine the deformation 
of the geometry caused by the fractional branes. 

It would also be interesting to determine whether the above computation coincides 
with calculating periods of the holomorphic three-form of the associated CY geometry. 
This problem will be left for future study. 



6. Conclusion 



We have seen in this paper that quantum deformed moduli spaces for geometrically en- 
gineered theories encode very non-trivial relations in the associated matrix models. The 
relations can be argued to be consistency conditions for the loop equations of the matrix 
model, when the loop equations degenerate. Indeed, the constraints are powerful enough 
to give a curve which solves the matrix model for a generic number of cuts. 

It is not clear under what circumstances the approach presented in this paper is valid. 
The most conservative possibility is that it applies only to those theories that correspond 
to a geometric realization of branes at singularities. This is suggested because the effects 
of the probe brane on fractional branes (what one would call the back reaction) cancels. 
It is very likely that this is tied to the fact that the low energy effective field theory of 
the probes on the moduli space is essentially AT = 4 SYM, and therefore it does not get 
corrected, nor does it develop a gaugino condensate. To this extent one can argue that their 
treatment should be essentially classical. This does not prevent the fractional branes from 
affecting the geometry of the probe, and we see that indeed they do by forcing geometric 
transitions. These transitions exactly deform the points in moduli space where the probe 
brane effective field theory has less super symmetry. 

The method serves to solve the theory for values of q which are roots of unity. Since 
in the SUSY gauge theory q is in principle a complex variable, one can expect that the 



-34- 



solution of vevs of chiral operators depends holomorphically on q, with possible singularities 
at some special values of q. These singularities should happen exactly when the associated 
SU(2) q representations constructed in the paper disappear (at roots of unity). For a given 
rank of the gauge theory TV, there are only classically a finite number of values of q which 
can be pathological (those that are r-roots of unity with r < N, or q = 0, oo). In principle 
this information should determine the pole structure of the holomorphic observables of the 
gauge theory as a function of q. Since we can in principle solve the theory for roots of unity, 
one should be able to use analyticity (holomorphy of the supersymmetric gauge theory) 
to solve the theory for any q, and address the issues of multi-cut solutions of the TV = 1* 
gauge theory. These and other related issues are currently under investigation. 
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Appendix 

In this appendix we treat the three matrix model of the paper by considering a saddle 
point for the eigenvalues of the matrix X. 
The potential of the theory is given by 



After one diagonalizes X, with eigenvalues Aj and if one integrates Y,Z, one obtains an 
effective potential for the eigenvalues of X given by 



This is effectively a one matrix model for X. Notice that in the limit where the 
eigenvalue distribution becomes continuous, for each cut in the eigenvalue plane associated 
to X, there are image cuts in the complex plane rotated by q and q^ 1 in the saddle point 
equation for the eigenvalues. For the SU (2)i representations one has eigenvalues on these 
images, and these give a picture of the SU(2) q representations. In the case of q a root of 
unity, there is a bound on the number of images per cut, and this makes the case of a root 
of unity special. 



7V/x _1 (tr (XYZ - qXZY) + tr (V(X))) 



i 



i<3 



The saddle point equations for the eigenvalues A are given by 




(6.2) 
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Also, one can derive the first quantum equation 3.46 for the case q n = 1 by multiplying 
the above saddle point equation by 

summing over all the eigenvalues, and dividing by N. However, the author has found no 
systematic way to extend this procedure to derive the other quantum equations from the 
above saddle point equations. 
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